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j^", We analyze the role played by entanglement in the dynamical evolution of Grover's search 

© : algorithm in the space of qubits. We show that the algorithm can be equivalent^ described 

as an iterative change of the entanglement between the qubits, which governs the evolution 
of the initial state towards the target state, and where the entanglement can be quantified in 
terms a single entanglement monotone. We also provide a necessary and sufficient condition 
for the quadratic speedup, which illustrates how the change in the bipartite entanglement 
of the state after each iteration determines the corresponding increase in the probability of 
i' finding the target state. This allows us to reestablish from the entanglement perspective 
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that Grover's search algorithm is the only optimal pure state search algorithm. 
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I. INTRODUCTION 

Entanglement — a consequence of the tensor-product structure (of subsystems) and superposition 
in quantum mechanics — has been slated to be one of the main protagonists in quantum computa- 
tion: since an efficient physical realization of a quantum computer necessarily requires it to have 

x 

' a scalable tensor-product Hilbert space — the space of qubits. 

a ■ 

Of course, entanglement is not necessary for quantum computation, since any computation in 
the space of qubits can be mapped onto a unary quantum system (for instance, an atom), where the 
computation will be devoid of entanglement — Hilbert spaces of same dimensions are fungible 
However, as clearly illustrated in Ref. [l], if a quantum computing system wants to avoid incurring 
exponential expense in the form of some physical resource for the computational problems which 
require exponential Hilbert space dimensions for its execution, then it is mandatory for the system 
to consist of subsystems (or the degrees of freedom) such that it's Hilbert space is equivalent to 
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the Hilbert space of qubits, i.e., the number of subsystems must grow linearly with the number of 
qubits required in an equivalent quantum computer. Nonetheless, the requirement of the scalable 
Hilbert space is necessary but not sufficient for a scalable quantum computer: it should further 
allow efficient implementation of the computational process itself — i.e., the initialization, control 
unitary dynamics, error corrections, and measurements. In a nutshell: entanglement leads to the 
saving of exponential spatial resources by accessing an arbitrary superposition of ./V states in just n 
(N = 2 n ) subsystems [2|, hence the motivation to understand it's temporal role (number of oracle 
queries required) in quantum computation. 

Jozsa and Linden [3] provided a major step towards understanding the role entanglement plays 
in speeding the dynamical evolution of a quantum computational process: a pure state quantum 
computation process necessarily requires multipartite entanglement which grows with problem size 
to achieve an exponential temporal speedup, and if the entanglement is capped to a fixed number 
of qubits — independent of the problem size — then the computation can be classically simulated 
with an equivalent amount of classical resources (also see this was demonstrated for Shor's 
efficient quantum factoring algorithm [fj]. 

This paper illustrates the temporal role of entanglement in Grover's search algorithm [6], and 
in the process further explains the result of Jozsa and Linden [3J]. Before we discuss our results, 
we first define the problem of searching a database. The problem is defined in terms of an oracle: 
one is given a very large unstructured database consisting of N (3> 1) items, one has to find a 
(multiple) marked item(s); the oracle consists of a function f(x): if one inputs an item x (an oracle 
query), it outputs 1 (yes) if x is a marked item, otherwise it outputs (no); to obtain a marked 
item with the minimum possible number of queries to the oracle is the search problem. In the space 
of qubits, Grover's algorithm [f], Q] executes the search by using only 0{^/N) oracle queries, but a 
classical digital computer — deterministic or probabilistic — would require 0{N) oracle queries on an 
average. Thus, quantum computers provide a quadratic temporal speedup over classical computers 
even though they both require 0(log 2 N) spatial resources to perform the search. What makes 
Grover's search algorithm significant is the fact that it is the only optimal [lj pure state search 
algorithm. 

Grover's search algorithm and the understanding of it which emerges from our paper can be suc- 
cinctly described as follows. The initial state \Sq), an equal superposition of N states representing 
an unsorted database with r marked states, is iteratively rotated towards the target state, where 
each iteration consists of the oracle operation followed by an application of the reflection operation 
(see Sec. Ull)- Each time the oracle operation inverts the marked states while leaving the unmarked 
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states unchanged, it creates a minimal structure in the state. This is exploited by the reflection 
operator by rotating the state about \Sq) to take one optimal step closer to the marked states. The 
structure is minimal in the sense that the dynamical evolution of initial state is restricted to the 
effective two dimensional space, i.e., the state after the kth iteration is given by 

\S k ) = A k \t) + B k \t ± ) , (1.1) 

where \t) represents an equal superposition of all the marked states, and \t±) represents the same 
for all the unmarked states. Thus, each iteration evolves the state simply by changing A k and B k , 
the two real parameters (excluding the normalization condition). 

The corresponding entanglement perspective of the algorithm (see Sec. IIIip : It is the iterative 
change of entanglement between the qubits that drives the initial state towards the target state. 
An iterative change arises because the oracle operation generates entanglement between all the 
qubits [l(3], which is then necessarily reduced by the corresponding reflection operator. The need 
for the iterative change in the entanglement can be motivated in the following way. The database 
consists of a large number of n-qubit states, and a small subset of which are the desired marked 
states — thus, each time the oracle generates the entanglement, it facilitates the corresponding re- 
flection operator to rotate the resulting state one step closer to the desired marked states. However, 
the consequence — or the limitation — of the dynamical evolution of the algorithm in the effective 
two dimensional space translates into the fact that the entanglement in the algorithm is restricted 
to a trivial bipartite form — i.e., it can be quantified by a single bipartite measure of entanglement. 
This is due to the fact that the Schmidt decompositions of the state after the kth iteration, \S)-), 
with respect to all the divisions of n qubits into two subsets, will consist of the Schmidt coeffi- 
cients which are necessarily functions of A k and B k (except for a constant factor). Now if you fix 
either A k or B k , it simultaneously fixes all the coefficients for all the decompositions. Moreover, 
since an entanglement measure is necessarily defined in terms of the Schmidt coefficients, an ar- 
bitrary Schmidt decomposition and a single measure of entanglement is thus sufficient to fix the 
entanglement of the state. 

More precisely, it is the change of entanglement in the state \Sk) due to the (k + l)th iteration 
that determines (Ak+i) 2 — (A k ) 2 , the change in the probability of finding the marked states from 
the kth to the k + 1th iteration. This follows from the equation derived in Sec. IIII CI 



where C(\Sk)) represents the concurrence 



111 ] — a measure of entanglement — of the state \Sk). We 



show that the above equation is a necessary and sufficient condition for the quadratic speedup, 
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and the integration of the equation, such that Ai changes from Aq to 1, determines the number 
of the oracle operations required for the search. This fact allows us to further reestablish from 
the entanglement perspective that Graver's search algorithm is the only optimal pure state search 
algorithm {J. 

The inference that can be drawn from our paper is that if a quantum algorithm requires 0(2 n / M ) 
oracle queries for it's execution, then it optimally exploits M effective dimensions of the 2 n (the 
problem size) dimensional Hilbert space — i.e., the pure initial state at all times during the evolution 
by the algorithm can be represented in terms of the same M orthogonal states. This will translate 
into the optimal exploitation of M-partite entanglement — i.e., M — 1 independent measures of 
entanglement will govern the evolution of the algorithm. Therefore, when one oracle query is 
needed for the execution of an algorithm (an exponential speedup), just as in Shor's algorithm [jj, 
then it would necessarily require optimal utilization of n-partite entanglement which grows with 
the problem size [3|. Simply put, the bipartite entanglement is the complete story of Graver's 
search algorithm — the reason for the quadratic versus the desired exponential speedup — which 
stems from the inherent inability of the oracle in generating any global structure. This implies a 
lack of multipartite entanglement between the qubits which grows with problem size, a necessary 
requirement for the exponential speedup [3]. In conclusion: the entanglement allows simultaneous 
saving of spatial and temporal resources when a quantum algorithm is executed in the space of 
qubits. 



II. GROVER'S ALGORITHM: THE SUPERPOSITION AND INTERFERENCE 

Here we give a brief summary of Grover's algorithm; for further details see We 
consider a database of ./V (3> 1) elements, and let it contain r (r < N) marked elements, which we 
want to find. The database is mapped onto the N states of a quantum system: 

\Xj) ; j = l,...N. (2.1) 

The first step of the algorithm is to form a equal superposition of the N states: 

1 N 



N . . 



Let us assume that r is known, and define the target state \t) as 



IT . 
J = l 
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a normalized linear combination of marked states. Similarly, the nontarget state \t±) represents 
the same for the unmarked states: 

1 N 

j=r+l 

By using equations (|2.3|) and (|2.4|) . |5q) can be reexpressed as 



= sin0|i) + cos 6"|t_L> 

= A \t) + B \t ± ) . (2.5) 

The state \Sq) is iteratively evolved to the target state \t), where each iteration consist of two 
unitary operations, the oracle operation Rq- 

Ro = 1- 2\t)(t\ , (2.6) 

followed by the reflection operator Rs : 

R So = 2\S )(S \-I . (2.7) 

By definition, the oracle Rq (|2.6p operation inverts the marked states, leaving the unmarked states 
unchanged; and the reflection operator Rs (|2.7p rotates the state about the hyperplane \So), hence 
the name 'reflection operator'. The central feature of the algorithm is that the iterative application 
of Rs Ro simply rotates \Sq) in the effective two dimensional hyper plane, {\t), \t±)}. This can be 
deduced by their action as described below. Let the state after kth iteration be 

\S k ) = A k \t) + B k \t ± ) . (2.8) 

An application of Ro on \S k ) gives 

\S k ) = -A k \t) + B k \t ± ) , (2.9) 
and a further application of Rs on the above state gives 



= A k+1 \t) + B k+1 \t ± ) . (2.10) 
By solving the recursion relation contained in Eq. (|2.10p we obtain 

A k = sin(2/c + 1)6* ; B k = cos(2/c + 1)6* , (2.11) 



6 



which can be verified by mathematical induction. Now, if we choose the number of iterations k 



as the nearest integer to (Tr/4)y/N/r, then the state \Sk) ~ \t), and a further measurement in 
the computational basis will provide one of the marked states. Moreover, the algorithm can be 
appropriately modified such that an arbitrary initial state can be used as the initial state, as 
long as (t\ift) = p is nonzero, to search the database in O(^i) oracle queries 13]. The obvious 
modification to the standard search algorithm (as described above) is to replace the reflection 
operator (|2.7p with the following reflection operator: i?u)/^,i = 2\ip)(ip\ — I. 

The search algorithm as presented in this section was independent of whether the database 
states (12. ip were of a quantum system, which may or may not be devoid of entanglement. 

Moreover, in principle, they could also represent the states of a classical system which allows 
superposition, for example, different modes of a classical electromagnetic wave. We refer the readers 
to Ref. [jjj] for the implementation of the oracle and the reflection operations in such systems; but, 
as discussed in the introduction, such an implementation necessarily incurs exponential spatial 
overhead. For our purpose here, we now map the algorithm onto the space of qubits. 

A. Grover's algorithm in the space of qubits 

The N database elements can be conveniently mapped onto the N = 2" product states of n 
qubits: 

\Xf) = \xi) ® . . . ® \x n ) ■ j = l,...N, (2.12) 

where \xi) G {|0), |1)} represents the states of the ith qubits — a two-dimensional Hilbert space Ti.2] 
and the superscript n denotes that \X™) G H® n . The database states (|2,12p are known 

as the computational basis. The initial state of the algorithm \Sq) (|2.2p is created by applying 
Hadamard transformation H [12] to each of the qubits in the product state |0)® n : 

(11 \ \ ® n N 

The oracle Ro (12. 6p is implemented via a conditional unitary transformation Uo on the computation 
basis states and an ancilla qubit \y), which is chosen to be in the state (|0) + |l))/2: 

Uo\X?)\y) = \Xp\y®f( Xj )) = (-l)fto)\X?)\y) , (2.14) 

where © denotes addition modulo 2, and f(xj) = 1 if \X™) is a marked state, otherwise f(xj) = 
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121 ] . The action of Uq from the computational basis perspective — ignoring the ancilla qubit from 
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the description, as it remains unchanged — reduces to Ro = I — 2\t n ){t n \ (|2.6p . The implementation 
of the reflection operator A2.TH follows from the Eq. (12.131) : 

R So = #®"(2|0 n )(0 n | - I)H® n = 2\S%){S%\ - I , (2.15) 

where |0 n ) = |0)® n ; and by construction it requires 0(n) gates. 

The description so far in this section implies that the search algorithm in the space of qubits 
requires 0{\J N/r) oracle queries to perform the search, and it's implementation requires 0{n) 
physical resources. In any situation where r is not known, one can either execute the algorithm 
by averaging over the several guesses for r, and still obtain a marked state in 0(w N/r) oracle 
queries; or first estimate r in time 0(a/ N/r) by a quantum algorithm provided in Ref. [3], and 
then execute the search. 



III. ENTANGLEMENT IN GROVER'S ALGORITHM 



Here we quantify the dynamical evolution of the entanglement between the qubits when the 
search algorithm is executed in the space of qubits. We illustrate how the change of entanglement 
generated after each iteration plays the central role in the search algorithm. To facilitate the 
discussion of the entanglement, and in order to quantify it, we first introduce a measure called the 
concurrence ll| ; the choice to use the concurrence was dictated by its simplicity (linearity) , and in 



any case, the essential results would remain the same under an another entanglement monotone [la ]. 



A. Concurrence 



We can always Schmidt decompose an arbitrary bipartite pure state |\E' ab ) of a Ttrja and TL D b 
dimensional subsystems: 

D 

|* a6 )=Ev^l*;>®l*i>, (3-1) 

where D = mm(D a ,D b ). The squared Schmidt coefficients fj,j are the eigenvalues of the reduced 
density operators, p a and p b , of the two systems, and the vectors and |^) make up the 

orthonormal bases that diagonalize the reduced density operators. If all but one of the Schmidt 
coefficients are zero, then the state is separable, otherwise the state is entangled. The concurrence of 
the bipartite pure state \ty ab ) is simply related to the purity of the marginal density operators 



C(\* ab }) = V2(l-tr[(p°)2]) 




(3.2) 



8 



C(\^ ab )) is conserved under local unitary transformations because it is a function only of the 
Schmidt coefficients. It varies smoothly from 0, for pure product states, to y/2(D — 1)/D, for 
maximally entangled pure states. 



B. The Schmidt decomposition of |S£) 

We want to quantify the entanglement of the state generated after the kth. iteration of the 
algorithm. Recall, the state after the kth iteration of the search algorithm is given by 

|5j?) = A k \t n ) + B k \tl) , (3.3) 

where (and henceforth) the superscript n denotes that the states are n-qubit states. We want to 
Schmidt decompose \Sj}) with respect to an arbitrary partition of n qubits into two subsets of (say 
the first) I qubits and the remaining (n — I) qubits; which will allow us to obtain the concurrence 
of the state by the use of Eq. (|3.2[) . 

For the sake of simplicity, we assume that there is a single target state (r = 1). We discuss 
the case of multiple marked states in the appendix IVH where we show that the results for a single 
marked state generalizes to the case when there are r marked states. To obtain the Schmidt 
coefficients corresponding to the bipartite decomposition of \S^), we define the target state to be 

\t n ) = \X[)\X?- 1 ) , (3.4) 

where \X\) is a I qubit product state (|2.12p . One can conveniently express the nontarget state as 

= — j = |V2^ - 1 1 x[\\N n - l \ + y/tf - 1 1 N l )\X?- 1 ) 

+ y/{2 1 - 1)(2™-* - l)|iV z )|JV n -')| , (3.5) 

where 

\N l ) = (2'-l)-5^|Xj) (3.6) 

2 n-l 

\N n - 1 ) = (2 n - 1 - 1)"3 \ X j~ 1 ^ ■ ( 3 - 7 ) 

3=2 

Equations (|3.4p and (|3.5p imply that the reduced density matrix pi obtained by tracing out (n—l)- 
qubit states from the density operator p 1 ^ = \S^){Sk\ can be represented in terms of the two 
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dimensional basis {|A{}, \N 1 }}, i.e., p l k (or l ) will have two non zero eigenvalues, which are the 
squared Schmidt coefficients of \S k ); and they are 

X l ± (k) = ^± \^~-r] 2 (A k -B k t a ne) 2 Blj , (3.8) 
where r/ is a constant that depends on the bipartite decomposition, 

r, = (N"- l \(N l \S%) = (^ (2 '~^ (2 _7~ 1) ) 2 • ( 3 - 9 ) 

and if N ^> 1, then r] ~ 1. Also, one can check that r\ is proportional to the entanglement of the 
nontarget state, i.e., C(|i™)) = 2r/ 

For the sake of completeness, we provide the Schmidt decomposition of \S k ): 



ffl = \J^k)W k ) + \r k - l u + J>L(k)W k )-w k ~ 1 )- , (3.10) 

where 

l4)± = *f\b\X[) + (A' ± (fc) - a)|JV')| , 

|$T'>± = AT'l^lXr') + (A±(A;) - a')| A™-') | , (3.11) 
where AA and A/"' are the normalization constants, and 

A 2 2 ™~' - 1 d 2 

a = A k + — — - B k , 
2' — 1 

« = A k + — — - B k , 



'2'-l (2"-'-l)V2^1 a 
^- T A fcJ B fe + _ B fc , 



2 n ~ l -l , (2 l - lW2 n ~ l - 1 o 

6 = V iv^r + " W^i Bk ■ (3 ' 12) 

The above Schmidt representation of |5^) makes an important point obvious, i.e., the main con- 
tribution to the Schmidt coefficients (|3.8p . and the corresponding Schmidt vectors, comes from the 
term A k B k . This arises from the entanglement between the target and nontarget states generated 
by the oracle — a fact central for the interpretation of our results, as discussed in the introduction. 
Now we use the Schmidt representation to analyze the dynamical evolution of the entanglement in 
the algorithm. 
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C. The concurrence of |S£) 

Equations (13.2|) and (13. 8p imply that the concurrence of \S k ) is given by 

C(|S2» = 2 v B k (A k -B k taa.e\ (3.13) 

= 2?7 sec sin 2fc0 cos(2£; + 1)0 (3.14) 

sec0 ( sin2fc0 \ dA\ 
~ V 29 ^sin(2A; + l)0 J dk { j 

- w • < 3 - 17 > 

where the first equality follows from (|3.2p and (|3.8p ; the second inequality is obtained by using 
A k = sin(2£; + l)0 and B& = cos(2/c + l)0 (12. lip . Moreover, the second equality assumes that we are 
in the first quadrant of the hyperplane {\t), \ t±)} (henceforth our discussion will be restricted to the 
first quadrant), which means that all the quantities in (|3.14p are positive; otherwise one has to put 
the absolute value sign in the right hand side of the (|3.14p , since by definition the concurrence (|3.2p 
is positive. The third equality (|A17|) is a different representation of the second, which follows from 
Eq. (|2.1ip ; the fourth equality (|3.16p follows from fact that sin 2k9 /sin(2/c + 1)0 ~ 1, and it is always 
positive but less than unity, except for the first few iterations; the fifth equality (|3.17p results from 
the assumption that N 3> 1 — the domain where the quadratic speedup is meaningful — then Aq ~ 9 
and rj sec 9 ~ 1, in which case the concurrence becomes independent of the bipartite decomposition. 

Equations (|3.13p - (|3.17p quantifies the entanglement between any I and n — I sets of qubits after 
the k iteration, except for a negligible constant factor. The Eq. (|3.17|) is the main result of the 
paper: it explains that the search algorithm exploits the change in entanglement after each iteration 
to evolve the initial state to the target state, i.e., the change in C(\S^}) due to the fe + lth iteration 
which determines (Ak+i) 2 — {A^) 2 , the change in the probability of finding the target state from 
kth to k + 1th iteration. Therefore, it follows that Eq. (|3.17p will not only governs the number 
of oracle queries needed to search the database, but, as we show in the next section IIVI that it is 
indeed a necessary and sufficient condition for the quadratic speedup. 

The evolution of the concurrence C(\SJi)) with the change in the number of the iterations k 
has been plotted in Fig 1: here we discuss a some of its salient features of it from the perspective 
of the change in the Schmidt coefficients with the change in k, further details can be obtained 
from Eqs. (|3TKj) - (|3T2l . First notice, Eq. (EHip implies CflSJ)) = (k = 0) and C(\t n )) = 
(k ~ \/iV(7r/4)), as it should be, since by choice |5q) and |i n ) are separable, otherwise C(\S^)) 
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is always nonzero. As k is increased, the difference between the Schmidt coefficients A+(A;) and 
X l _(k) starts decreasing, therefore C(\S^)) monotonically increases, and attains it's maximum value 
when k m ^/N(tt/8), where X l + (k) = XL(k) = l/y/2. Further iterations monotonically decreases 
C(\S%}), since the difference between the Schmidt coefficients again starts to increase. When 
C(\Sfc)) approaches zero, it signals that the target state is being approached, i.e., when X + (k) = 0, 
X l _(k) = 1 (or A k = 1, B fc = 0). 

We now give an alternate description of the change in C(\S^)) with a change in k, which is 
more relevant for the purpose of this paper. This description naturally arises from the entangling 
properties of the oracle and reflection operator. 




1000 



FIG. 1: Plots of C(\S%)) against k for r = 100 and N = 10 s (r? = 1). 



D. Entangling and Disentangling process in the algorithm 

We first show that the oracle operation of selectively inverting the target state generates the 
bipartite entanglement. This can be seen from the concurrence of the state Ro\S)-)- 

C(R \SD) = 27 1 (A k + B k t a n6)B k (3.18) 

= 2r ? sec6»sin2(A; + 1)0 cos(2A: + 1)6* (3.19) 
_ s_ec9 ( sin2(fc + l)fl \ cL4| 

" V 26 [sm(2k + l)8 J dk ^ Uj 

where the second equality (|3.19|) is obtained by substituting the explicit values of A k and B k given in 
pTTT]) : the third equality (f33U1) also follows from Eq. ([2lT]) . where sin2(A: + l)6»/sin(2A; + 1)6 » 1, 
except for the first few iterations, but because sin2(A; + l)0/sin(2/c + 1)6 > sin 2k6/s'm(2k + 1)6, 
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then Eqs. (lA17|) and (ET201) imply that d(R \S%)) > C^S^)). This can be also seen directly from 

C 1 (i? |5^))-C 1 (|5^)) = 4r ? i?ftan# (3.21) 

= 4i] cos 2 (2k + 1)6 tan 2 9 , (3.22) 

which (in the first quadrant) is always positive but a decreasing function (see Fig. 2). 

The corresponding reflection operator increases the probability of finding the target state, then it 
should reduce the entanglement between the qubits. This can be seen by evaluating the difference, 
Ci(|5^ +1 )) — Ci(Ro\S%)). The difference is always negative and decreasing, as shown in Figure 
2 (where we have plotted the negative of the difference). Fig. 2 provides a simple explanation 
for Fig. 1: the monotonic increase of C(\Sk)) as K is increased from k = to k = w N/r(ir/8) 
results because, in this range of k, each oracle operation generates more entanglement than the 
reduction due to the corresponding reflection operator. Whereas, between k = sj N/r(n/8) and 
k ps v / iV/V(7r/4) the converse occurs. More importantly, we now show that the iterative increase 
and decrease in entanglement due to the oracle and reflection operator, respectively, is optimal. 
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FIG. 2: Plots of C(R \S k )) - C(\S k )) (solid line) and C(R \S k )) - C(R So Ro\S k )) (dotted line) against k 
for r = 100 and N = 10 8 . 

IV. QUADRATIC SPEEDUP: THE NECESSARY AND SUFFICIENT CONDITION 

If an unsorted database of N items has a single marked item, then Grover's search algorithm 
finds the marked item in 0(y / N/r) oracle queries. Then, we have already shown, that the state 
after the kth iteration (or oracle query) \S%) has entanglement given by Eq. (|3. 17p ; therefore, it 
is a necessary condition to search the database. To show the converse, assume that there is some 
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(arbitrary) search algorithm which employs the initial state \Sq) (|2. 13|) to search the database, and 
it searches in the hyper plane then our previous discussion implies that the Schmidt 

decomposition of the n-qubit pure state produced after each iteration with respect to an arbitrary 
division of n qubits into two subsets will involve only two nonzero Schmidt coefficients, excluding the 
normalization condition. Let the state generated by the arbitrary algorithm after the kth. iteration 
be 15^)', and let = sin 2 (f>(k) and \^ = cos 2 4>(k) be its squared Schmidt coefficients. Now, 
assume that the concurrence of \Sk}' is given by (l/2Ao)d(Ak) 2 /dk, where A\ is the probability of 
finding the target state in \S^)', and Aq = sin# = sm(p(k = 0) ~ <j). This implies 



1 d{Ak? - C(\S k 



n\/\ 



2tj) dk 



= sm2(f)(k) , (4.1) 

where the second equality follows from the definition of the concurrence ()3.2f) . Integrating the 
above equation, and imposing the initial condition Aq = <fi, implies A\ = sin 2 (2/c + 1)0, hence 
the condition C(\S%)') = (l/2Ao)d(Ak) 2 /dk is necessary and sufficient to achieve the quadratic 
speedup. 

A natural question presents itself: If an algorithm searches out of the two dimensional hyper- 
plane, can it do better than than the quadratic speedup? The search executed out of the hyperplane 
can happen if the oracle instead of inverting the target state introduces a relative phase between 

nn 

the target and nontarget state — see for example, Grover's fixed point algorithm [16. 1171] . Then it 
is easy to show that it will imply C(\Sk}) < (l/29)d(Ak) 2 /dk, because the projection in the hyper 
plane will always reduce the concurrence [la ], therefore it will provide no temporal advantage over 
Grover's search algorithm. Moreover, one can apply an operation which is a more general opera- 
tion than the reflection operation, i.e., it is not restricted to the effective two dimensions of the 

I 

search algorithm. Then, as shown explicitly by Zalka [21], nothing additional can be gained: the 
shortest distance between the target and nontarget state, dynamically speaking, lies in the two di- 
mensional hyperplane {\t), i.e., along the geodesic — hence the reflection operator is optimal. 
This can also be argued from the entangled perspective as follows: although the entanglement of 
the state which rotates in a multidimensional hyperplane would involve more than one measure of 
entanglement — it won't be fixed just by the concurrence — but it is the oracle operation that guides 
the evolution to the desired target state, and by definition, it is optimally restricted in the two 
dimensional hyperplane. 

To show that, as far as the pure state search is concerned, Grover's search algorithm is the 
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only optimal search algorithm, we have to show that the quadratic speedup is also asymptotically 
optimal g], i.e., for any number of oracle queries, and for arbitrary reflection operators. 



A. Asymptotic optimality 

Our proof of the asymptotic optimality is analogous to the proof in Ref. 8]. Suppose we waul, 
to search an unsorted database of N items with a single marked item. We map the problem onto a 
N dimensional space of n (2 n = N) qubits. Consider an arbitrary initial n-qubit pure state \ip n ), 
which contains a single marked state \t n ). \ip n ) is evolved by invoking the oracle, R l = 2\t)(t\ — I, 
as follows: 

\<j%) t = UtR'oUt-iRo ■ ■ ■ U^oW 1 ) . (4.2) 

where the unitary operators Uj's are arbitrary, and the probability of finding the target state after 
T oracle queries is given by Alp 2 = | (t n \(f) T )t\ 2 ■ The idea of the proof is to compare the above 
evolution to the case when |^ n ) is evolved without invoking the oracle Rq: 

\4%) = U T ...U 1 \i> n ) , (4.3) 

in which case, let the probability of finding the marked state after T oracle queries is given by 
At 2 = |(i n |<^j.)| 2 . The main part of the proof works by obtaining an upper bound for the difference 
\Aip 2 — At 2 \, and averaging over N linearly independent choices of the marked states, since one 
can always design a special algorithm which is suited for a particular choice of the marked state. 
The upper bound can be obtained by considering Eq. (|3.17|) : 

< >/2, (4.5) 
where Aq 2 is the probability of finding the marked state in the initial state \ip n ), and the second 



inequality is obtained from the fact that the concurrence is bounded by yj2{N — 1)/N w y/2 (see 

secims]). 

Integration of the above equation (|4.5p provides the upper bound: 

J2\Al 2 - A k 2 \ < V2TVN, (4.6) 
t 

where we have used Aq = y/l/N. Now to complete the proof one obtains the lower bound by 
considering a worst case scenario: we want to invoke the oracle large enough times such that we 
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should be able to distinguish sufficiently via a measurement all for N linearly independent 



moreover, this should be true for all N choices of the marked states: 

N 

J2\A{ 2 - A k 2 \ > eN . (4.7) 
t=i 

Then, gSJ) and (j47)) implies T > O(VN). 

Therefore, quadratic speedup is also asymptotically optimal |8|]. Thus, Grover search algorithm 
is the only optimal algorithm for searching an unsorted database with a pure state. The question: 
Can parallel quantum computation improve the quadratic speedup? 

V. ENTANGLED PARALLEL QUANTUM COMPUTING 

Here we show that when Grover's search algorithm is executed by I entangled computers — or 
/ different sets of n qubits — then it provides a certain linear advantage, i.e., they can produce I 
copies of the target state, where each copy has r marked states, with just 0(y / N/r) oracle queries. 
The advantage can be motivated in the following way. The result of the search algorithm is an 
equal superposition of all the marked states, and a further measurement in the computational basis 
provides one of the marked states. Suppose each of r marked states encodes different information, 
therefore, to extract all the information, the search algorithm needs to be executed (say) I times 
to produce I copies of the target state, which will require 0(ly/ N/r) oracle queries. 

We now show that the search with a multipartite entangled state can provide the linear ad- 
vantage. We denote the computational basis states of the feth quantum computer as 
j = 1 . . . N, and the initial state of the I quantum computers is a generalized GHZ state: 



choices of \t n ), which in turn implies that for some large k, and for some fixed e, A l k 





3 



(5.1) 



Now define the target state of I computers as 




(5.2) 



and similarly, the nontarget state \ can be defined: 




(5.3) 
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By using equations (|5.2p and (|5.3p . \Sq)i can be rewritten, just as in the standard search algorithm 
(see Sec. HI|) : 

= sin6\t}i + cos 6 i\t±) 

= Aol^i + Bolt^t . (5.4) 

Now, the standard search algorithm is executed by a single computer (say, the first), and the rest 
of of computers do nothing, i.e., the oracle operation R l for the I entangled computers is made 
of the product of I unitary operations, where the first unitary operation is the standard oracle 
operation (|2.6|) . and the rest of the unitaries are the respective identity operators. Similarly, the 
reflection operator R 1 Sq for the computers can be defined. It is straight forward to show that after 
m = (tt/4) yN/r iterations the local density operators of all the computers reduces to 

;El*?WI- ^ 

Notice, that the linear advantage is achieved via the entanglement in the initial state, although the 
joint unitary operations are local. 

VI. CONCLUSION 

Our discussion in this paper was restricted to the search algorithm, where the initial pure 
state was an equal superposition of all database states; however, it can be easily generalized to an 
arbitrary pure initial state. Moreover, if the search is performed with a mixed initial state, then 
our result, Eq. (|3.17p . has to be optimized not only over all the ensemble decomposition of the 
state, but also over all reflection operators. This, we will consider in a future publication. 

Since the entanglement in the search algorithm is of the limited bipartite form, it explains why 
the quadratic (temporal) speedup can be achieved if the algorithm is implemented in a classical 
system which allows superposition, as explicitly shown by Lloyd [14j. This is further reflected by 
the fact that even a classical digital computer can trivially simulate the algorithm. However, as the 
number of bipartite entanglement measures required for a quantum computational process grows 
(the multipartite- ness of entanglement increases), then the temporal expense incurred by the simu- 
lation of the computational process by a classical system, irrespective of whether the superposition 
is allowed or not, will necessarily grow; more importantly, the simulation will have no efficient 
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realistic description [1] of the computational process in the limit where the multipartite entangle- 
ment needed for the computations grows with the problem size fS]. After all, what distinguishes 
quantum computation — whether there is entanglement or not — from classical computation — where 
superposition is allowed or not — is that an arbitrary quantum state and it's dynamics has no effi- 
cient realistic description: the reason, perhaps, that makes quantum computation inherently more 
powerful than classical computation [1(. 

The main contribution of this paper is an explicit illustration of the dynamical role of en- 
tanglement in the search algorithm: the change of entanglement after each iteration governs the 
evolution of the initial state. Although the entanglement in the algorithm is limited, it is opti- 
mally exploited by the search algorithm. It is precisely due to the limited bipartite entanglement 
in the algorithm, that the algorithm leads to the quadratic speedup. The bipartite entanglement 
itself is a consequence of the oracle: which by definition, optimally restricts the evolution of the 
algorithm to the two effective Hilbert space dimensions. The simplicity of the algorithm does 
succeed in illuminating that entanglement — perhaps, entangling operations in the case of a mixed 
state computation — plays an essential role in saving both the spatial and temporal resources when 
exponential Hilbert space dimensions are required for the execution of a quantum computational 
process. Thus, entanglement is indispensable for a scalable quantum computer. 



APPENDIX A: MULTIPLE MARKED STATES 



Here we show that all our previous analysis of the role of entanglement in Grover's search 
algorithm, where we had assumed that there is a single marked state, can be generalized when 
there are r arbitrary marked states, i.e., the target state is given 

n = 4l>?>> ( A1 ) 

where \X™) are computational basis states (|2. 12j) . This implies that \SJ^) may not be separable in 
the arbitrary bipartite decomposition of the n qubits. Recall, the state after the fcth iteration is 
given by 

\St) = A k \t n ) + B k \tl) , (A2) 



and we want to obtain it's concurrence corresponding to the bipartite division of n qubits into I 
qubits and n — I qubits. Let's assume, without any loss of generality, that the r marked states 
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belong in the set of database states: 

M={|Xj)|X£-')} ; j = l,...,p, k = l,...,q, (A3) 

where p < r, q < r, and pq > r; and |Aj) and \X™~ 1 ) are /-qubit and (n — /)-qubit computational 
basis states, respectively. To obtain the concurrence of \S^}, we first define a convenient basis: 

i T ') = ^Ei x i)' ( A4 ) 
i^> = -^=Ei^j>. ( A6 ) 

v 2 - P j=i 

IN"'*) = 7 =_ g IX-') , (A7) 

where |T') and |T n ~') are a sum of all |Xj)'s and \X"j~ )'s, respectively, which are contained in M; 
similarly, |JV') and |A" n_ ') are sum of all LYj)'s and |X"~')'s, respectively, that are not in M. 
Notice, we can always express \t n ) as follows: 

3=1 

where Y^jPj = r > an d 

i p j> = ^X>j>> ( A9 ) 

where the sum is over pj number of states |Aj) which belongs in M, and {P l j\P l k ) = Pjk/ \/PjPk- 
The nontarget state by the use of Eqs. (|A4j) - (|A7|) . can be conveniently expressed as 

1*1) = -^=(^|r')|r^)-^|r) 



+ ^q(2 l - p)\N l )\T n ~ l ) + ^{2 l - p)(2»-' - q)\N l )\N n ~ l ) j . (A10) 
By substituting (jSS]) and fATOj) in (|A2j) . one obtains 

\Sft) = (A k - B k t a ne)j2\l^\Pt)\Xr l ) 



» r 

3=1 



+ -J±= Up-q\T l )\T«~ l ) + ^(2< - p)\N l )\T n - 1 ) 



+ ^{2l-p)(W-l-q)\N l )\N n - l ) \ , (All) 
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where we have substituted tan# = y / r/(N — r). By tracing out the (n — /)-qubit states from the 
density operator p£ = \S^)(S^\, one can obtain p l k ; which can then be substituted in the expression 



C7(|fl£» = «/2[l-tr (pi) 2 ], (A12) 



to show that the square concurrence of \S£) is given by 

C\\SD) = ^'^'^ (A, - B k t^6fBl (A13) 
+ (Al-Blt a n6) 2 C 2 (\t n )), (A14) 

where the first term represents the entanglement generated by the search algorithm after the fcthe 
iteration, but the second term is simply a byproduct of the initial entanglement, since it is not 
affected by the oracle operation. Therefore, when analyzing the role of entanglement in the search 
algorithm, it is necessary to ignore the second term in the above equation. If we do so, then the 
concurrence of the state \S£) generated by the algorithm is given by 

C(\S%)) = 2r,'(A k - B k t<m6)B k , (A15) 

where 

rf = {N- l \{N l \St) = ( (2 '~^"~'~ g) ) 2 , (A16) 

and when there is a single target state, i.e., when C(\t n )) = 0, p = q = 1, and r( — > r], it reduces 
(so does Eq. (|A14p to Eq. (|3.13p . Thus, the dynamical evolution of the algorithm, when there are 
r marked states, in terms of the concurrence can be similarly expressed as 

C(\S n )) - n' — ( S[n2k ° (All) 
C[lbk)) ~ V 26 ^sin(2fc + l)0 ) dk (A17) 

1 dAl . 

(A18) 



2A Q dk ' 

where Aq = sin 6 = y / r/A r . Below, we illustrate the derivation of Eq. (|A14|) via a simple example, 
where there is a little loss of generality, and it is more instructive. 

1. An entangled target state 

Suppose we Schmidt decompose \S%) with respect to the division of n qubits into n — 1 qubits 
and the remaining nth qubit. There is no loss of generality here, because it is sufficient to quantify 
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the entanglement of \Su), otherwise it would mean that more than two parameters and (ex- 
cluding the normalization condition) is being changed after each iteration. Moreover, the Schmidt 
decomposition of |5^) with respect to (n — l)-qubit nth-qubit states can be used to prove that 
Eq. (|A18P is necessary and sufficient for the quadratic speedup when there are multiple marked 
states, just as we did for a single marked state (see Sec. II Vj) . 

Now, suppose that the target state \t n ) in the bipartite decomposition is given by 

1 



where 



, Vp\P n ~ l ) ® |0) + ^\Q n ~ l ) ® |1> ) , (A19) 



IP™ -1 ) = — V IX"- 1 ) , (A20) 



io B - l > = ^=Ei-xr 1 >, (A21) 

where p + q = r. We further assume for the sake of simplicity that \P n ~ l ) and \Q n ) are 
orthogonal, i.e., Eq. (|A19P is the Schmidt decomposition of \t n ). If |P n_1 ) and \Q n ~ 1 ) are not 
orthogonal, which will generally be the case in the computation basis, then one can start with the 
decomposition of \t n ) as given in Eq. (|A9p . \t n ) ()A19|) can be conveniently reexpressed as 

\t n ) = sin0 \P n ~ l ) tg> |0) +COS0 IQ™- 1 ) ® |1) , (A22) 

where sin0 = \fpfr and cos</> = -J q/r are the Schmidt coefficients. If either p or q is zero (or 
= 0), then the target state is separable, otherwise it is entangled. The entanglement in \t n ) is 

C(|i n )) =sin20 , (A23) 

obtained by using Eq. (|3.2p . This implies that \t r \) can be decomposed as 
1 



= - 7 =^|P«-i)|l) + ^|Q-i)|0) + yiV72^7|iV«-i)(|0) + |l))j 

= tan^^sin^lP™- 1 )!!) +cos0|Q"- 1 )|O)^ + \A - tan 2 ^| A^ 1 ) ( ^jp ) j , (A24) 
where |A^™ _1 ) is defined to be 

JV/2 

==£iAr- 1 ), ( A25 ) 



I TV"- 1 ) 



v r+l 
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therefore, by construction, \N n 1 ) is orthogonal to \P n 1 ) and \Q n 1 ). The concurrence of the 
nontarget state is given by 



C(\tl)) = tan 6^2 -tan 2 6(1 + cos 2 20) . (A26) 
The initial state \Sq) (|2.13p is by definition separable in all the bipartite decomposition: 

\SS) = \ V2 sinO sin ^\P n ~ l ) + v^sinflcos^Q"-') + V^OlN 11 - 1 }^ , (A27) 



which is obtained by substituting ()A22p and (|A24|) in Eq. (fA~2]h 

The state after fcth iteration, by using equations (|A22p and (|A24|) . can be expressed as 

\S£) = AJ S in0|P^)|O) + cos0|Q^)|l) 

+ B k ttm6\Bim<f>\P n - l )\l) + cos (/>|Q n ^)|0) J 

+ B k Vl-t a ^e^U n - l ) m ^ ll)) ^j . (A28) 

The two eigenvalues of the reduced density operator of the nth qubit, obtained by tracing out the 
states \P n ~ l ), \Q n ~ l ) and \U n ~ l ), can be used in Eq. (|3.2p to show that the square concurrence of 
| S%) given by 

Cl(\St)) = Cl(\S n k )) + C 2 2 (\S n k )) , (A29) 

where 



C x {\Sl)) = 2(A fe -B fc tan0)£ fc J(l-tan 2 0) (A30) 



= 2(A k -B k tan8)B k ^ v - _ - ' , (A31) 

and 

C 2 {\SD) = (A 2 -5 2 tan#)sin20 (A32) 

= (A 2 k - B 2 k tan 9)C(\t n )) . (A33) 

Notice, Ci(\S k )) and C2(\S k )) has the same form as in Eq. <|A14|) . Ci(\S k }) is term which is the 
concurrence generated by the search algorithm, and has been discussed in Sec. lIII Al We restrict our 
attention to the properties of C^dSJ?)); it is zero for the initial state (the initial state is separable by 
definition), otherwise nonzero everywhere else; it monotonically increases and attains its maximum 
value of sin 2(f) — the initial concurrence of the target state. Thus, C^dS^)) is simply a byproduct of 
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the initial entanglement of the target state. The above example give us the opportunity to illustrate 
the special, but interesting, case r = N/4 ([7]). It highlights in a rather more dramatic fashion the 
entangling and disentangling nature of the oracle and the refection operator, respectively. 

a. r=N/4 

When r = N/4, then it requires just one oracle query to search the database — classically, it 
would require two oracle queries. If r = N/4, then Eq. ([23]) implies Aq = 1/2 and B = y/3/2, 
thus the initial state expressed in the bipartite form in Eq. (1A27|) reduces to 

\So) = -j=(a*<t>\P) + cos0|Q) + \U)^j r^^-J (A34) 

Then the action of Ro on the above state generates a maximally bipartite entangled n — 1-qubit 
and nth qubit state: 

Rom = i=^in0|P^)-cos0|Q"-^ + ^^(^TT 1 ) • ( A35 ) 

i.e., Ci(Ro\Sq)) = 1. Again notice, it is the the oracle which creates entanglement between the 
target and nontarget states. Now the action of Rs on the above state reduces it to the target 
state: 

R So Ro\S%) = \t n ) , (A36) 

then C\ (Rs Ro \ So)) = 0; thus, the reflection operator reduces entanglement in the search algo- 
rithm. 

ACKNOWLEDGEMENTS 

The author thanks Apoorva D Patel for useful discussions on Grover's search algorithm, and 
quantum computation in general. 



[1] Robin Blume-Kohout, Carlton M. Caves, and Ivan H. Deutsch, Journal of Optics B: Quantum and 

Semiclassical Optics 6, S801-S806 (2004). 
[2] A. Ekert and R. Jozsa, Phil. Trans. Roy. Soc. Lond. A, 356, 1769 (1998) 
[3] Richard Jozsa and Noah Linden, Proc. R. Soc. Lond. A, 459, 2011 (2003) 



23 



G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 

P. W. Shor, in Proceedings of the 35th Annual Symposium on Foundations of Computer Science, edited 
by S. Goldwasser (IEEE Computer Society, Los Alamitos, CA, 1994), pages 124-164; P. W. Shor, SIAM 
Journal of Computing, 26, 1484 (1997). 
L. K. Grovcr, Phys. Rev. Lett. 79, 325 (1997). 

M. Boyer, G. Brassard, P. Hoyer, and A. Tapp, Fortsch. der Phys.-Prog. of Phys. 46, 493 (1998). 

C. H. Bennett, E. Bernstein, G. Brassard, and U. V. Vazirani, SIAM Journal of Computing, 26, 1510 

(1997). 

C. Zalka, Phys. Rev. A 60, 2746 (1999). 

S. L. Braunstein and A. K. Pati, Quantum Information and Computation, 2 399 (2002). 
P. Rungta, V. Buzek, Carlton M. Caves, M. Hillcry, G. J. Milburn, Phys. Rev. A 64, 042315 (2001). 
M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge Uni- 
versity Press, Cambridge, England, 2000). 
L. K. Grovcr, Phys.Rev.Lctt. 80, 4329 (1998) 
S. Lloyd, Phys. Rev. A 61, 010301 (2000). 
Guifrc Vidal, J. Mod. Opt. 47, 355 (2000). 
L. K. Grovcr, Phys. Rev. Lett. 95, 150501 (2005). 

T. Tulsi, L. K. Grovcr and A. Patcl, Quant. Inform, and Comput. 95, 150501 (2006). 
Pranaw. Rungta (to be published). 



